Dopant metrology in advanced FinFETs 
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Ultra-scaled FinFET transistors bear unique fingerprint-like device-to-device differences attributed to random single 
impurities. This paper describes how, through correlation of experimental data with multimillion atom tight-binding 
simulations using the NEMO 3-D code, it is possible to identify the impurity's chemical species and determine their 
concentration, local electric field and depth below the Si/Si02 interface. The ability to model the excited states rather than 
just the ground state is the critical component of the analysis and allows the demonstration of a new approach to atomistic 

impurity metrology. 

PACS numbers: 



I. RECENT PROGRESS IN DONOR 
SPECTROSCOPY 

Modern transistors are getting so small that it is in- 
creasingly difficult to use traditional techniques for their 
study and their characterisation [T]. This is particu- 
larly true for the identification of the impurity chemical 
species present in the channel and for the quantification 
of their concentration. Different groups have recently in- 
vestigated the effects of ultra-scaling in silicon Field Ef- 
fect Transistor (FET) geometries and interesting results 
have emerged: 

• The fBM Zurich group [Ij has shown that screen- 
ing due to interface traps in ultra-scaled silicon 
nanowires can cause substantial increase on the ion- 
ization energy of the dopant. This result has pro- 
found implications for the design of future FET de- 
vices. 

• M. Pierre et al [3] and R. Wacquez et al |4j have 
shown that the presence of a single dopant in the 
channel of a frigate FET can dramatically alter its 
electrical signature, even at room temperature. 

• M. Fueschsle et al |5_, have studied the band struc- 
ture effects on single-crystal silicon geometries for 
which the source, the drain and the gates are fabri- 
cated using an atomically sharp doping procedure 
[6]. Therefore they have investigated on the conse- 
quences of scaling a device up to the single atom 
limit. 

• M. Tabe et al 7J have demonstrated that, in ultra 
scaled silicon FET devices, even in the presence of 
a dopant rich environment, it is possible to observe 
the signature of a single dopant. 
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Overall, all these studies have allowed a better un- 
derstanding of the effects that arise due to the ultra- 
scaled environment, however, they have also indicated 
that, for the successful design of future Complementary- 
Metal-Oxide-Semiconductor (CMOS) devices, a substan- 
tial amount of knowledge is still missing. As an example 
it is not yet clear how CMOS technology will be able to 
overcome critical challenges such as scaling-induced vari- 
ability of device characteristics [H [S] . 
As a consequence, it is of interest to discuss more in de- 
tails a method that can be used to demonstrate atomic 
impurity metrology 110). In fact, through correlation 
of experimental data with multimillion atom simulations 
in NEMO 3-D, the impurity's chemical species can be 
identified and their concentration, local electric field and 
depth below the Si/Si02 interface can be determined 
[H [To]. Furthermore, the extension of the dopants in 
the source/drain (S/D) regions can be measured by spec- 
troscopy of confined states in the channel [HI [TU] . Lastly, 
the effective current distribution in the channel can be 
determined by means of thermionic emission theory [llj . 
Following these lines, the goal of this paper is to intro- 
duce a new dopant metrology technique to be used for 
ultra-scaled CMOS devices. 



II. TRANSPORT BASED DOPANT 
METROLOGY IN ADVANCED FINFETS 

Direct or impurity-mediated tunnelling between source 
and drain competes with the thermally activated current 
and thereby affects the sub-threshold swing. Advanced Si 
MOSFET devices have shrunk to dimensions where the 
magnitude of the sub-threshold swing is dominated by 
the nature of individual impurities in the channel region. 
Due to the random nature of the impurity distribution, 
parameters like excited band levels [5] , ionization energy, 
threshold voltage and leakage current show device-to- 
device fluctuations [2H11 13 111 fT^HTi] . An atomistic point 
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FIG. 1: Geometry and electrical characteristics of a single donor located in the channel of a FinFET device, (a) Scanning 
Electron Micrograph of a typical FinFET device, (b) Band diagram along the x-direction with the D'^-state in resonance 
combined with the measured source/drain current versus gate voltage for a typical sample. QD^ and QD^ indicate resonances 
of a quantum dot, formed by the confinement provided by the corner effect and residual barriers in the access regions between 
source/drain and channel. The gate voltage where the band edge in the channel is aligned with the Fermi energy Ep in 
source/drain, indicated by Ecb, is estimated by subtracting one unit of addition energy from QD^. Below the band edge, there 
are resonances ascribed to the D" and D~ charge states of a single donor. 



of view is imperative to understand and determine the 
underlying donor characteristics in the channel. Several 
recent experiments have showed that the fingerprint of a 
single dopant can be identified in low-temperature trans- 
port through such devices [31 21 [7| jTUl [T5HT7] , suggesting 
a new method to characterise devices down to the level of 
a single impurity. In the context of quantum state control 
in quantum computing, it was shown that it is possible 
to model the orbital levels of a single donor in the chan- 
nel of scaled FinFETs by means of multi-million atom 
modelling [lOl [T8j . In this chapter it is shown that the 
aforementioned technique can be used to perform single 
donor mapping in advanced Si FinFETs. By carefully ex- 
amining the fingerprints of isolated donors in an ensemble 
of devices, the chemical species, concentration, and local 
field of donors in the channel region can be determined 
in a non-intrusive fashion. 



III. DEVICES 

The FinFET devices used in this study consist of 
crystalline silicon wires (fins) with large contacts pat- 
terned by 193 nm optical-lithography and dry etching 
from Silicon-On-Insulator, see Fig. [l] (a). After a boron 
channel implantation, a 100 nm polycrystalline silicon 
layer was deposited on top of a nitrided oxide (1.4 nm 



equivalent Si02 oxide thickness). A phosphorus (P) im- 
plant was used for pre-doping and the structure was sub- 
sequently patterned using an oxide hard mask to form 
a narrow gate. Next, high-angle arsenic (As) implanta- 
tions were used for source or drain extensions, while the 
channel was protected by the gate and 50 nm wide ni- 
tride spacers and remained p type. Finally, As and P im- 
plants and a NiSi metallic silicide were used to complete 
the source or drain electrodes. The samples described in 
this chapter all have a gate length of 60 nm. 

Transport measurements are performed on an ensem- 
ble of devices at a temperature of 4 K and a search for the 
fingerprints of isolated donors is performed. These single 
donors are located in or near the active cross-section of 
the channel, i.e. the cross-section of the FET body where 
the potential is lowest and the electrical transport thus 
takes place [TT] . Large electric fields induced by the gate 
or even corner effects can reduce the active cross-section 
to dimensions much smaller than the FET body |11| . 
Furthermore, corner effects play a major role in these 
devices, as indicated by an active cross-section of only 4 
nm^, determined by thermionic transport measurements 

A single donor's fingerprint is characterised by a pair of 
resonances in the source-drain current, Isd, versus gate 
voltage, Vq, characteristics at low Vsd (Fig- [l] (b)). The 
positions of a pair of resonances in Vg are an indication 
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FIG. 2: Source/drain difFerential conductance of the D'' charge state as a function of bias voltage and gate voltage of typical 
single donor FinFET devices. The excited states, indicated by the black dashed lines, form the fingerprint by which the donor 
properties can be identified. The dots are a direct indication for the energy of these states {Ei = eVso , with Ei the i-th excited 
state and e the unit charge.) a) Sample 13G14: Excited states are observed at 3.5, 15.5 and 26.4 meV b) Sample 10G16: 
Excited states are observed at 2, 15 an 23 meV. c) Sample GLJ17: Excited states are observed at 2, 7.7 and 15.5 meV. 



of the energy of the one-electron (D°) and two electron 
charge states (D~). A large quantum dot present in the 
channel (with charge states indicated by QD^ and QD^) 
is also observed, which allows a rough determination of 
the position of the band edge in the active area. While 
the quantum dot in the channel is almost always found, 
only about one out of seven devices shows the fingerprint 
of a donors. The identification of the resonances of the 
donor is based on the determined binding energy, charg- 



ing energy and the odd-even spin filling |17| . Next, the 
excited energy levels of the one-electron (D°)-state are 
determined by sweeping both the Vg and Vsd biases and 
measuring the differential conductance {dlsD/dVsD) in 
the appropriate bias space, see Fig. [2] In this so-called 
stability diagram the typical diamond-shaped region as- 
sociated with Coulomb-blocked transport between the 
and D~ states is observed. The total electronic transport 
in the conducting regions increases as an excited level of 
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FIG. 3: a) Calculated wave-function density of an As donor with d = 4.3 nm for three different fields. The gray plane indicates 
the Si02-interface. From low-fields (where the donor has a bulk-like spectrum) to high fields, the donor electron makes a 
transition from being localised on the donor to being localised at the silicon interface, b) First eight eigen-levels of an As donor 
3.2 and 4.3 nm below the interface and a P donor 3.2 nm below the Si02 interface as a function of electric field (F) calculated 
in a tight-binding model (NEMO 3-D). Note that we measure excited states relative to ground state (lower black line in this 
graph.) 



the D^-state enters the bias window defined by 
source/drain, giving the stability diagram its character- 
istic pattern [11] [T7] indicated by the dashed black lines. 
The red dots indicate the combinations of Vsd and Vg 
where the ground state is at the Fermi energy of the 
drain and an excited state is at the Fermi energy of the 
source. It is the bias voltage Vsd in this combination 
that is a direct measure for the energy of the excited 
state (eVsD.N = En), where En is the energy relative to 
the ground state and N is the level index). The excited 
states as determined in this fashion are depicted in Table 
I. These levels are not bulk-like but heavily influenced by 
the local electric field and the nearby interface [18ti20j . 



Finally, the measured level spectrum is compared to a 
multimillion atom tight-binding NEMO 3-D calculation 
pi] of the system taking two possible chemical species. 
As and P, into account. NEMO 3-D solves for the eigen- 
values of the single electron Schrodinger all band equa- 
tion in a tight-binding approach [22] . The NEMO simu- 
lation package is based on about 14 years of development 
at Texas Instruments, NASA Jet Propulsion Laboratory, 
and Purdue University [211 El] . Each atom is explicitly 
represented and the electronic structure of the valence 
electrons is represented by ten sp3d5s* orbitals. Spin 
can be included explicitly into the basis by doubling the 



number of orbitals. Spins are coupled through spin-orbit 
coupling resulting in accurate valence band states. The 
five d orbitals help shape the curvature of the conduction 
bands to achieve appropriate masses and symmetries at X 
and L. The tight-binding parameters are tuned to repro- 
duce the bulk silicon properties under various strain con- 
ditions faithfully. For systems where the primary interest 
is in the conduction band properties and if no magnetic 
fields need to be considered, spin can safely be ignored 
without any significant loss of accuracy. Effects due to 
crystal symmetry, strain, local disorder, and interfaces 
can be explicitly included in the model through direct 
atomic representation. 

The single impurity states are modelled with a sim- 
ple Coulomb potential away from the impurity site and 
a central on-site core correction to match experimen- 
tally observed bulk-like impurity energies. The simu- 
lation domain for a bulk-like single impurity must be 
large enough such that the hard wall boundary condi- 
tions imposed by the finite simulation domain are not 
felt by the central impurity. With a simulation domain 
of 30. 4x30. 4x30. 4nm'' corresponding to about 1.4 million 
silicon atoms the impurity eigenstates move less than 1 
/xeV with further domain size increases. A critical mod- 
elling capability here is the need to be able to compute 
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Device 




El (meV) 


E2(meV) 


E3(meV) 


d (nm) 


F (MV/m) 


s 

(mcV) 


10G16 


Exp 


2 


15 


23 






TB 


2.2 


15.6 


23.0 


3.3 


37.3 


0.59 


11G14 


Exp 


4.5 


13.5 


25 






TB 


4.5 


13.5 


25.0 


3.5 


31.6 


0.04 


13G14 


Exp 


3.5 


15.5 


26.4 
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26.3 


3.2 


35.4 


0.17 


HSJ18 


Exp 
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10 


21.5 






TB 


4.5 


9.9 


21.8 


4.1 


26.1 


0.63 


GLG14 


Exp 


1.3 


10 


13.2 






TB 


1.3 


10 


12.4 


5.2 


23.1 


0.28 


GLJ17 


Exp 


2 


7.7 


15.5 






TB 


1.3 


7.7 


15.8 


4.9 


21.9 


0.77 



TABLE I: First three measured excited states of each sample (see for example Fig. 2) versus the best fit to the NEMO 3-D 
model (as depicted in Fig. 3). The fit yields a unique combination of (F, d) for each single donor device. The measurement 
error for each level is estimated to be around 0.5 meV. 



reliably the ground states as well as the excited states of 
the impurity system, as that is a significant component of 
the impurity fingerprint. Fig. |3] (b) shows typical eigen- 
state spectra for an As donor (two donor depths) and 
a P donor as a function of electric field. Three electric 
field regimes can be distinguished (Fig. |3] (a)). At the 
low field limit (F mV/m) the spectrum of a bulk As 
donor is obtained. In the high field limit (F ~ 40 MV/m) 
the electron is pulled into a triangular well formed at 
the interface and the donor is ionized [TOl [HI HO]- In 
the cross-over regime (F ~ 20 MV/m) the electron is 
de-localised over the donor- and triangular well poten- 
tial and the level spectrum consists of levels associated 
with the donors, levels associated with the triangular well 
at the gate interface (formed by the local field) and hy- 
bridised combinations of the two (TUl HSl HH] ■ 

The measured level spectra are least-square fitted to 
a calculation over a sufficiently sized region of F-d (field 
and donor depth parameter space). At least three ex- 
cited levels per donor are taken into account to make 
the fit over-determined. The fitting procedure is per- 
formed for two different species of donor atoms. As and 
P, which were both used in the fabrication process. The 
concentration of donors in or near the active area of the 
FET channel can be estimated by comparing how many 
times a single donor is identified with the relevant volume 
where donors can be found. 



IV. RESULTS 

Of the forty-two devices that have been examined, six 
have been found to exhibit the fingerprint of a single 
donor in the transport characteristics. These devices 
were subsequently measured carefully and the D° level 



spectrum was fitted. The quality of the fits, indicated by 
X2, across the six samples is 0.92 (x2 <1 means a good 
fit) assuming As donors. For P donors, we find a x2 of 
10.22 (although two of the six are comparable in qual- 
ity.) Based on the fits, the donors active in the devices 
are assumed to be As. The (over-determined) fit fur- 
thermore yields a unique combination of (F, d) for each 
single donor device, as shown in Table I, together with 
the measured level spectra, and their fits. 

Finally, it is also possible to obtain the donor concen- 
tration. As mentioned before, the active cross-section of 
the devices are heavily reduced by a corner effect |11) . 
The identified donors are either in or near these corners, 
and from the least-square fit we actually know the donor 
depths range from ~ 3 to 6 nm below the Si/Si02 in- 
terface, see Table I. The part of the FinFET channel 
where donors are found is a volume spanned by the donor 
depths and the gate length (60 nm), and a single donor 
device is found in average only in one out of seven de- 
vices that are measured. This means that for each ~ 1000 
nm^ of measured material it is possible to found one As 
dopant, and, as a consequence, a local As concentration 
of about lO^^cm"-^ can be extrapolated. 



V. CONCLUSIONS 

In this paper a novel method for dopant metrology 
is introduced. The excellent quantitative agreement be- 
tween the measured and modelled level spectra gives an 
indication of the high level of confidence of this method 
to determine the chemical species and local field of single 
impurities in silicon FinFET transistors. Furthermore, 
the local concentration of donors near the active cross- 
section of the FinFET can be estimated. The present 
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method offers opportunities for non-invasive characteri- 
sation down to the level of a single donor and could be 
a future tool in the guidance of device processing. This 
is especially true for such CMOS devices for which vari- 
ability problems are increasing dramatically with ultra- 
scaling and this justifies the importance of the method 
described here. 
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